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C^ : Abstract 

We consider the m-phase Whitham's averaging method and pro- 
pose the procedure of " averaging" of non-local Hamiltonian structures. 
The procedure is based on the existence of a sufficient number of local 
commuting integrals of the system and gives the Poisson bracket of 
^ I Ferapontov type for Whitham's system. The method can be consid- 

ered as the generalization of the Dubrovin-Novikov procedure for the 
, , local field-theoretical brackets. 

>< 
^ ■ Introduction. 

We consider the averaging of the non-local Hamiltonian structures in 

Whitham's averaging method. As well known Whitham's method permits 

to obtain the equations on the "slow" modulated parameters of the exact 

periodic or quasi-periodic solutions of systems of partial differential equations 

and it was pointed out by Whitham ([!]) that these equations inherit the local 

Lagrangian structure if the initial system has it. The Lagrangian formalism 

for Whitham's system is given in this approach by the "averaging" of the 

local Lagrangian function for the initial system on the corresponding space 

of (quasi)-periodic solutions. Some basic questions concerning Whitham's 

method can be found in [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. 



B.A. Dubrovin and S.P. Novikov investigated also the question of the 
conservation of local field-theoretical Hamiltonian structures in Whitham's 
method and suggested the procedure of "averaging" of local field-theoretical 
Poisson bracket to obtain the Poisson bracket of Hydrodynamic type for 
Whitham's system ([4, 7, 9], see also [13]). 

The Jacobi identity for the averaged bracket and the invariance of the 
Dubrovin- Novikov procedure of averaging was proved by the author in [14] 
(see also [15]) using the Dirac restriction procedure of the initial bracket 
on the subspace of quasi-periodic "m-phase" solutions of the initial system. 
The connection between the procedure of Dubrovin and Novikov and the 
procedure of averaging of the Lagrangian function in the case when the initial 
local Hamiltonian structure just follows from the local Lagrangian one can 
be found in [16]. 

Some extension of the averaging of " local" Hamiltonian structures for the 
case of discrete systems is also present in [17]. 

In the present work we deal with the Poisson brackets having the non-local 
part of the form 



fc>0 

+ Y. (^kS\k){v, Vx, ■ ■ ^J^^x - y)Sf,^-^{<^, ^y, . . .) (0.1) 

A:>0 

where Ck = ±1, i^{x — y) = —I'iy — a;), dx^^x — y) = 5{x — y) and both 
sums contain the finite numbers of terms depending on the finite numbers of 
derivatives of ip with respect to x. 

We point out here that the brackets (0.1) can be found in the so-called 
"integrable" systems (see [18, 19, 20]). 

The most general form of the non-local Hamiltonian operators (0.1) con- 
taining only 5'{X — Y) and 5{X — Y) in the local part and the quasi-linear 
fluxes S)^-.AU)Ux of " hydrodynamic" type in the non-local one was suggested 
by E.V.Ferapontov in [23] as the generalization of the bracket introduced in 
[22]. We shall discuss here the possibility of "averaging" of the brackets (0.1) 
in Whitham's method to obtain the bracket of such "Hydrodynamic type" 
for Whitham's system. 

As was shown by Ferapontov the Hamiltonian operators of this type re- 
veal a beautiful differential-geometrical structure following from the Jacobi 



identity of the bracket ( [23, 24, 25, 26]). (In particular they can be obtained 
as the Dirac restriction of local differential-geometrical Poisson brackets on 
the space with fiat normal connection ( [24]).) 

The first example of the non-local bracket (of Mokhov-Ferapontov type, 
see [22]) for Whitham's system for NS equation in the one-phase case was 
constructed by M.V.Pavlov in [34] from a nice differential-geometrical con- 
sideration. After that there was a question about the possibility of construct- 
ing of the nonlocal Hamiltonian structures for Whitham's system from the 
structures (0.1) for the initial one. As was mentioned above the Hamiltonian 
operators (0.1) exist for many "integrable" systems like KdV and in the work 
[26] (see also [35]) there was a discussion of the possibility of averaging of the 
non-local operators for KdV equation using the local bi-Hamiltonian struc- 
ture and the recursion operator for the two averaged local Poisson brackets. 
The corresponding calculations for the m-phase periodic solutions of KdV 
were made by V.L. Alekseev in [36] . 

Here we propose the general construction for the averaging of operators 
(0.1) in Whitham's method which is the generalization of the Dubrovin- 
Novikov procedure for the case of the presence of non-local terms in the 
bracket. Our procedure does not require the local bi-Hamiltonian struc- 
ture and can be used in the general situation. Like as in the procedure of 
Dubrovin and Novikov we require here the existence of the sufficient number 
of commuting local integrals generating the local flows according to (0.1) and 
we also impose the conditions of "regularity" of the full family of m-phase 
solutions as in the local case (see [14]). 



1 Some general properties of the non-local 
brackets. 

Let us consider the non-local 1-dimensional Hamiltonian structure of type: 
{V\x), V^{y)} = Y: i?;i)(^, ^., . . ■)S^'\x - y) + 

fc>0 

+ J2^{k)iv,Vx,---Mx-y)flf^){<^,<^y,...) l<i,j<n (1.1) 



fe>0 



where we have the finite numbers of terms in the both sums depending on 
the finite numbers of derivatives of ip with respect to x. 

We shall call the local translation invariant Hamiltonian function the 
functional of the form: 



H[if]= JVHiv,^.,-.-)dx (1.2) 

Here z/(a; — y) is the skew-symmetric function such that 

D^v{x - y) = 6{x - y) (1.3) 

and 5^^^ {x — y) is the /c-th derivative of the delta-function with respect to x. 
We assume here that the bracket (1.1) is written in the "irreducible" 
form, that is the number of terms in the second sum is the minimal possible 
and the sets S'(^^ and TL\ are both linearly independent. From the skew- 
symmetry of the bracket (1.1) it follows then that the sets of S)j^^ and Tl-. 
coincide and it can be easily seen that the bracket (1.1) can be represented 
in the "canonical" form 



{^\x), ^^{y)} = Y: B^v, v., ■ ■ ■)S^'\x - y) + 

k>0 

+ H efc5(\)(v9, v?^, . . :)u{x - y)Sli^^{^, <py, . . .) (1.4) 

A:>0 

where Cfc = ±1 . 

Indeed since the sets {5*/^)} and {T^} coincide we have the one finite- 
dimensional linear space generated by fiuxes (vector fields) 

and the symmetric (view the skew-symmetry of the bracket and the function 
z/(a; — y)) finite-dimensional constant 2-form which describes their couplings 
in the non-local part of (1.1). So we can write it in the canonical form 

according to its signature after some linear transformation of the fiows S(^k) 
with constant coefficients. 

We should also define in every case the functional space where we consider 
the action of the Hamiltonian operator (1.4) and this can depend on the 
concrete situation. The most natural thing is to consider the functional space 



(p{x) and the algebra of functionals I[ip] such that their variational derivatives 
multiplied by the flows S(k){'^, fx, ■ ■ •) give us the rapidly decreasing functions 
as \x\ ^ oo. Here we shall use the functionals of type 

/n 
J2 <^^{x)qp{x)dx 
p=i 

where qp{x) are arbitrary smooth functions with compact supports to get 
all the properties of the bracket (1.4) and for the other functionals used in 
the considerations we shall assume that they have the compatible with the 
bracket (1.4) form in the sense discussed above. 



We construct here the procedure which gives us the brackets of Ferapon- 
tov type ( [23]- [26]) 

{WiX), U^{Y)} = g''^{U)6'{X -Y)+ h'{^{U)U^5{X - Y) + 
+ E e,Sl,^^{U)U^^v{X - Y)S^,^,iU)U^ I < u, ^i, \,6 < N (1.5) 

A;>0 

from the initial bracket (1.4) of the general form after the averaging on the 
appropriate family of exact m-phase solutions of local Hamiltonian with re- 
spect to (1.4) system with the local Hamiltonian function H. 

So for the Poisson brackets (1.4) we consider Whitham's method for the 
local fluxes generated by Hamiltonian functions (1.2) (if they exist), i.e. 

^l = Q\^,^,,...) (1.6) 

Certainly this can be only for special Hamiltonian functions in general 
case so all the considerations here appeal as a rule for the "integrable sys- 
tems" like KdV, NS, etc. where we have a lot of such functionals. 

We now formulate some general theorem about the non-local part of the 
bracket (1.4). 



Theorem 1.1 

Suppose we have the non-local Hamiltonian operator written in the "canon- 
ical" form (1-4), 



Then: 

1) The flows 

0^ = ^J,)(^,^.,...) (1.7) 

commute with each other. 

2) Any of the flows (1.7) conserves the Hamiltonian structure (1.1) 

Proof. 

Let us consider the functional 






:ii 



for some qp{x) with the compact supports and consider the Hamiltonian flow 
^*(a;) generated by (1.8) according to (1.4), i.e. 



k>0 '^•^ ^ A:>0 



X 



SLi^P'^Pz, ■ ■ .)qj{z)dz - / SIA^,^^, . . .)qj{z)dz 



'(ky 



Cfc = ±1 (there is a summation over the repeated indices). 
For the Hamiltonian flow ^^{x) we should have: 



:i.9) 



C.J 



{x,y) = 



:i.io) 



where J is the Hamiltonian operator (1.4) and C^ is the Lie-derivative given 
by the expression: 



C,J 



'\x,y) = Jc{^)j^^JHx,y)dz- 



J''{z,y)-4r^e{x)dz- fj^^{x,z)-4--e{y)dz 



6(p'^{z) J 5(p'^{z) 

Let us now consider the relation (1.10) for x and y larger than any z from 
the support of qp{z). Then we shall have 



[C^jf (x, y) = Y: bUv, v., ■ ■ ■)S^'Hx -y) + 



k>0 



+ H ekSkJ<^, v^x, • • -Mx - y)SL){'^, i^y, . . .) + 



fe>0 



+ H efe5(\)(v9, <^x, ■ ■ -Mx - y)Sl^){v, fy, ■ ■ •)- 



A:>0 






X 



if 



Slki){'^,Vw,---)qp{w)dw 



fc>0 "-^ \A:'>0 



if 



Slk'){v,'^w, ■ ■ ■)qp{w)dw 

•' k>0 

6 / ■ r 1 /■°° 

E ^fe'r^TT -^(V)!^' ^x, • • •) - / Sf^J^, v?^, . . .)qp{w)dw 
J dzJ2 ekSlk)i(p, (Px, ■ ■ -Mx - z)S-lk){ip, ipz, ...)x 



k'>0 



k>0 



X 






1 /"OO 

Z J—oo 



where B]^ {f, ^Px, ■ ■ •) and SLJip, (px, ■ ■ •) are the derivatives of these functions 
with respect to the flow 



Vl 



Y.^kS\k){V,'P 



k>0 



Sfk){<p,'^ni,...)qp{w)dw 



:i.ii) 



Here we used that x,y > Supp qp in the expression (1.9) for ^{x) and 
^■^{y) and also when omitted the variational derivatives with respect to (p^{x) 
and ip''^{y) of the non-local expressions containing the convolutions with qp{w). 
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So we have 

()=[C^j\\x,y) = 



k>0 



k>0 



Y, efc5(fc)(v^, Vx, ■ ■ ■)i^ix - y)Sl){<^, <^y, . . .) 



fe>0 






X 



if 

Z J— oo 



'5ffe/)(v5>V^«;,---)?p(^)c^^ 



fc>0 






1 r"^ 
2J-C 



k>0 



A;'>0 



'^'5(V)(v^>v^^>---) 



5v9*(2;) 



Z J— oo 



- dzJ2 efc5(fc)(v^, <^x, • • O'^la^ - ^)'5ffc)(v5, V^^, . . .)x 

•^ fc>0 



E fife 



fc'>0 



6ip^{z) 



Z, J— CO 



'5(l./)(v^,^«„---)?p(^^)c?^i^ 



- dzJ2 ekSlk)i^, (fz, . . -Mx - y)Sl^^{ip, cpy, ...)x 

•' fc>0 



Y (^k'Slk>){<^,<^x,---)7^ 



k'>0 



J dzJ2 (^kSlk)i(p, (Px, ■ ■ -Mx - y)Slk){ip, (pz, . . .) 



)qp(w)dw 



X 



.+ 



Yl (^k'SLJ^,^y,...) 



k'>0 



dif'^ 



^^j_jUv.V^....)qA 



w)dw 



k>0 



1 r°° 

2 J ^^"^W 



Sfi,)i(p, (p^a,, ■ ■ ■)qpiw)dw 



CkJ 



{x,y) + 



■ J2 ekek'Slf,,)i(p,(px,---)Sff^-^{ip,ipy,...)x 

k,k'>0 



1 



'(fc')l 

6 



^2i(^(^)(^'^---->V(.) 



6 



S'^^^.,^{(p,(p^,...)qp{w)dw 



5(V)(v9, Vz, ■ ■ O r ,/ X [J_^ S'^(k)iv, Vw, ■ ■ .)qp{w)dw 



)dz 



;i.i2) 



where 

(1.7) 



ChJ 



V 



X, y) are the Lie derivatives of J with respect to the flows 



As can be easily seen the last term in (1.12) is the only one containing the 
functions non-equal to zero when x ^ y which are not skew-symmetric as a; — >■ 
y (the other non-local terms contain the function z/(x — y) ). So from (1.12) 
we have that it should be identically zero for any qp{w) with the support 
satisfying the conditions described above {x,y > Supp qp{w) ). Using the 
standard expression for the variational derivative and the integration by parts 
we obtain that this term can be written in the form 



Y. ekek'Slk,)i(p, (f^, . . ■)SlJip, (py,... 



k,k'>0 



qpi^) 



S{k), S(^k') 



[z)dz 



for any qp{z) {x,y > Supp qp{z)) where [5'(fe), S^k')] is the commutator of the 
flows (1.7). So for the linearly independent set of Si^k) we obtain 



'S'(fc), 5'(fc/) 







From (1.12) we then have also for the linearly independent set of Si^k) and 
different qp{w) that 



9 



CkJ {x,y) = 

So we obtain the statements of the theorem. 

Theorem 1.1 is proved. 

It is also obvious that the statements of the theorem are vahd for all the 
brackets (1.1) written in the "irreducible" form since all S(k) and T(^k) in this 
case are just the linear combinations of the flows S(^k) ■ 



Remark. 

Let us point here that the first statement of the Theorem for the non-local 
brackets (1.5) of Ferapontov type was proved previously by E.V. Ferapontov 
in [23] using the differential-geometrical considerations. In [23]- [26] also 
the full classification of the brackets (1.5) from the differential geometrical 
point of view can be found. 

It is easy to see now that the local functional of type (1.2) 

V{ip,ipx,-- ■)dx 

generates the local fiow in the Hamiltonian structure (1.1) if and only if the 
derivative of its density V{(p,ipx, ■ ■ •) with respect to any of the fiows (1.7) 
is the total derivative with respect to x, i.e. there exist such Q(k){'^, '^x, ■ ■ •) 
that 

Vr^{^,Vx,- ■ •) = dxQ{k){V,Vx,- ■ •) 

or, as was pointed by E.V.Ferapontov ( [23]), this means that the integral / 
is the conservation law for any of the systems (1.7). 

From the Theorem 1.1 we have now that the flows (1.7) commute with 
all the local Hamiltonian fluxes generated by the local functionals (1.2) since 
they conserve both the Hamiltonian structure and the corresponding Hamil- 
tonian functions. 
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2 The Whitham method and the "regularity" 
conditions. 

Now we come to Whitham's averaging procedure (see [1]- [10]). Let us 
remind that in the ?72-phase Whitham's method for the systems (1.6) we 
make a rescahng transformation X = ex, T = et to obtain the system 

ev9^ = g*(v9,e^x,eVxx,...) (2.1) 

and then try to find the functions 

S{X,T) = {S\X,T),...,S"^{X,T)) 
and 

oo 

<^\e,X,T,e) = J2e'^<!>l,){e,X,T) 

k=0 

- 27r-periodic with respect to each 6" {9 = (0^, . . . , 9^)) such that the func- 
tions 

^\9, X, T,e) = f^ e'^^\,^{9 + ^^^^^, X, T) (2.2) 

fc=0 ^ 

satisfy the system (2.1) at any 9 in any order of e. 

It follows that $(Q-) {9, X, T) at any X and T gives the exact m-phase 
solution of (1.6) depending on some parameters U = {U^, . . . , U^) and initial 
phases 9q = {9l, . . . , 9^) and, besides that: 

S^ = tu"(f/(X, T)) , SI = A;"(f/(X, T)) 

where u^lU) and k°'{U) are respectively the frequencies and the wave num- 
bers of the corresponding m-phase solution of (1.6). 

The conditions of the compatibility of system (2.1) in the first order of e 
together with 

give us Whitham's system of equations on the parameters U{X, T) 
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U!f = y;(f/) U^ (2.3) 

which is the quasi-hnear system of hydrodynamic type. 

The first procedure of averaging of the local field-theoretical Poisson 
brackets was proposed in [4]- [9] by B.A.Dubrovin and S.P.Novikov. This 
procedure permits to obtain the local Poisson brackets of Hydrodynamic 
type: 

{[/^(X), U^{Y)} = g''^'{U)5\X -Y)+ V^^'iU) Ul 5{X - Y) (2.4) 

for Whitham's system (2.3) from the local Hamiltonian structure 

k>0 

for the initial system (1.6). 

The method of Dubrovin and Novikov is based on the presence of A^ 
(equal to the number of parameters U^ of the family of m-phase solutions of 
(1.6)) local integrals 



r = J V''iip,ip^,,...)dx (2.5) 

commuting with the Hamiltonian function (1.2) and with each other 

{F,H} = , {r,/^} = (2.6) 

and can be written in the following form: 

We calculate the pairwise Poisson brackets of the densities V^ in the form 

{V''ix),V^iy)} = Y: An^, cp,, . . .)<5W(x - y) 

k>0 

where 

A7iv, Vx, ■■■) = d^Q^^i^, ^^, . . .) 

according to (2.6). Then the Dubrovin-Novikov bracket on the space of 
functions U{X) can be written in the form 

{W^iX), U^{Y)} = {A7){U) 6'{X -Y) + ^-^Ul 6{X - Y) (2.7) 

12 



where (. . .) means the averaging on the family of m-phase solutions of (1.6) 
given by the formula: ^ 

1 f^ 

and we choose the parameters Lf^ such that they coincide with the values of 
I" on the corresponding solutions 

f/- = (P^ix)) 

The Jacobi identity for the averaged bracket (2.7) in the general case was 
proved in [14] (for the systems having also local Lagrangian formalism there 
was a proof in [16]). 

Let us note also that we get here the Poisson bracket only if we average 
the initial Hamiltonian structure on the full family of m-phase solutions (see 
[13], [14]). More precise requirements we shall formulate when describe the 
averaging procedure in the non-local case. 

Brackets (2.4) can be described from the differential-geometrical point of 
view. Namely, for the non-degenerated tensor g''^^ we have that it should be 
a fiat contravariant metric and the values 

should be the Levi-Civita connection for the metric g^fi (with lower indices). 
The brackets (2.4) with the degenerated tensor are more complicated but 
also have a nice geometrical structure (see [12]). 

The non-local Poisson brackets (1.5) are the generalization of local Pois- 
son brackets of Dubrovin and Novikov and are closely connected with the 
integrability of systems of hydrodynamic type reducible to the diagonal form 
( [21]). Namely, the system reducible to the diagonal form and Hamilto- 
nian with respect to the bracket (1.5) satisfies (see [23]- [26]) the so-called 

"'^ Strictly speaking this formula is valid for the generic set of the wave numbers /c", 
but we should use in any case the second part of it for the averaged quantities to obtain 
the right procedure. (Here /c" are continuous parameters on the family of the ?7i-phase 
solutions). 
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" semi-Hamiltonian" property introduced by S.P.Tsarev ( [21]) and can be in- 
tegrated by Tsarev's "generalized hodograph method". In [29] can be found 
also the investigation of the equivalence of "semi-Hamiltonian" properties 
introduced by Tsarev and the Hamiltonian properties with respect to the 
bracket (1.5). 

Let us also point here that the questions of integrability of Hamiltonian 
with respect to (1.5) systems which can not be written in the diagonal form 
were studied in [30]- [33]. 

The procedure of averaging of the non-local Poisson brackets in Whitham's 
method and the proof of the Jacobi identity for the averaged non-local bracket 
resemble the same things for the local brackets. However the difference in 
formulas of averaging and in the proof contain very many essential things and 
also many considerations valuable for the local case can not be used in the 
non-local one. So we should make here the consideration for the non-local 



case. 














The 


m 


-phase 


solutions 


of (1.6) 












V'' 


'ix,t) = 


^'{ujt + kx + 9o 








where 


UJ ■■ 


= {u^\.. 


.,u;"^),k = {k\ 



are the 27r-periodic solutions of the system 



uj^'^go. - Q\<l>, r$ec, r A;'^<l>e«e/3, . . .) = , a = l,...,m (2.9) 

depending on uj and k as on the parameters. So we assume that we have 
from (2.9) for the generic u and k the finite dimensional subspace Aico,k (in 
the space of 27r-periodic with respect to each 6*" functions) parameterized 
by the initial phase shifts 9q and maybe also by some additional parameters 

iy-4 J- iy-4 ft' ^ 



^For the multiphase case (m > 2) it is essential that the closure of any orbit generated 
by the vectors (wi, . . . , Wm) and (k\, . . . , km) in the 0-space is the full m-dimensional torus 
T™. For the case of "rationally-dependent" wi, . . . , ujm and ki, . . . ,km and to > 3 we have 
that the operators (2.9) are independent on each of such closed submanifolds in T™ which 
can have the dimensionality < m. The functions from Mu;.k can be found in this case 
from the additional requirement that they also the ?7i-phase solutions for (1.7) (with some 
Q?^-,(w, fc,r)) and the systems generated by the functionals I'^ (see later) also with some 
uj°"^{ijj,k^r). All these requirements uniquely define the finite-dimensional spaces AA^j^k 
which continuously depend on the parameters llj and k. 

14 



Combining all such Ai^j^k at different u and k we obtain that the m-phase 
solutions of the system (1.6) can be parameterized by A^ = 2m + h param- 
eters U^, . . . , U^ - invariant with respect to the initial shifts of ^" and the 
initial phase shifts 6^ after the choice of some "initial" functions $(j„)(6', U) 
corresponding to the zero initial phases. The joint of the submanifolds M.ui,k 
at all uj and k gives us the submanifold M. in the space of 27r-periodic with 
respect to each 9'^ functions which corresponds to the full family of m-phase 
solutions of (1.6). 

For Whitham's procedure we should now require some "regularity" prop- 
erties of the system of constraints (2.9). Namely 

(I) We require that the linearized system (2.9) at any "point" of Ai^^k 
has exactly h + m = N — m solutions ("right eigen vectors") C,(q)uj,k{(^,r) for 
the generic uo and k given by the vectors tangential to M.u),k which are the 
derivatives ^ga[6, r, u, k) and ^ri{d^ f^, ^, k) (at the fixed values of uo and k). 

(II) We also require that the number of the linearly independent "left 
eigen vectors" /t(g)a;,fc(^,T) orthogonal to the image of the introduced linear 
operator is exactly the same (N — m) as the number of the "right eigen 
vectors" ^{q)to,k{6,i^) for the generic oo and k. In addition, we shall assume 
that H(q)uj^k{(^,i^) o/so depend continuously on the parameters U" on M.. 

The requirements (I) and (II) are actually connected also with Whitham's 
procedure and the asymptotic solutions (2.2). Indeed, it is not very difficult 
to see that every fc-th term in the expansion (2.2) is determined by the 
defined above linear system with the nontrivial right-hand part depending 
on the previous terms of (2.2). We suppose that this system is resolvable 
on the space of 27r-periodic with respect to each ^" functions if the right- 
hand part is orthogonal to all the "right eigen vectors" /t(g)^fc(^,r) for the 
corresponding {oj = St, k = Sx,r). If so, the solution of this system can be 
found modulo the "left eigen vectors" C,{q)ui,k{(^,f^) with the same u, k and r 
defined by the zero term of (2.2). The corresponding freedom can be used 
then in higher orders of e to get the compatibility for higher orders systems 
defining the terms in (2.2) (see [2]). The compatibility conditions in the first 
order of e together with 

kT = OOx 
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give us the Whitham's system of equations (2.3). ^ 

We now discuss the requirements (I) and (II) from the Hamiltonian point 
of view. 

First of all for the procedure of averaging of the bracket (1.4) we need the 
set of the integrals f^, z/ = 1, . . . ,N like in the procedure of Dubrovin and 
Novikov, i.e. satisfying the following requirements: 

(A) Every I" is the local functional 



r = jv^{^,^^,...)dx (2.10) 

which generates the local flow 

v\. = Q\,){v,Vx,...) (2.11) 

with respect to the bracket (1.4). 

As was pointed above we should have for this that the local flows (1.7) 
deflned from the bracket (1.4) in the "canonical" (or "irreducible") form con- 
serve all the f^, i.e. the time derivatives of the corresponding Vl^p, ^Px, ■ ■ •) 
with respect to each of the flows (1.7) are the total derivatives with respect 
to X 

^V{v, Vx, . . .) = dxF};,^{v, vx, ■ ■ .) (2.12) 

for some functions Fn^Jip, ipx, • • •)• 

(B) All /^ commute with each other and with the Hamiltonian function 

(1.2) 

{r,F} = o , {r,H} = o (2.13) 

(C) The averaged densities (V) 



■^We assume that $(fc)(6',X, T) depend continuously on X and T and so are also 
well-defined for the non-generic fci, . . . , km and wi, . . . , ujm in the multi-phase Whitham's 
method. 
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can be regarded as the independent coordinates U^, . . . , U^ on the family of 
m-phase solutions of (1.6)."^ 

From the requirements above we immediately obtain that the flows (2.11) 
commute with our initial system (1.6) and with each other. 

From the Theorem 1.1 we obtain also that the commutative flows (1.7) 
defined by the Poisson bracket also commute with (1.6) and (2.11) since 
they conserve the corresponding Hamiltonian functions and the Hamiltonian 
structure (1.1). 

Now we can consider the functionals 

r = lim^ r V''iip,ip^,...)dx (2.15) 

c^oo 2c J-c 



and 

H=\im^j-^£vH{v,Vcc,---)dx (2.16) 

on the space of the quasiperiodic functions (with m quasiperiods). 

It is easy to see now that the local fluxes (1-6), (1-7) and (2.11) being 
considered on the space of the quasiperiodic functions also conserve the val- 
ues of I^ and H and commute with each other since these properties can he 
expressed just as the local relations containing (/?, ipx, ■ ■ ■ and the time deriva- 
tives of the densities Vl^p, ipx, ■ ■ ■), Vni'^, Vx, • • •) (it the same point x. 

Now we can conclude that all the fluxes (1.7) and (2.11) leave invariant 
the family of m-phase solutions given by (2.9) and can generate on it only 
the linear shifts of the initial phases 6q which follows from the commutativity 
of the flows 

VU{9) = S(\)(v9,A;>e.,A;"A;^v^,ce,,...) (2.17) 

and 

Vl^ie) = Ql^{^,k^^e'^,k^k('^e^g,, . . .) (2.18) 

with the flows (fla = (p^a and 

^Here again we should use everywhere the second part of the formula (2.14) for the 
averaged values on A4 for the right procedure as will be shown later. 
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on the space of 27r-periodic with respect to each 9°^ functions and the conser- 
vation of the functionals J*^ (i.e. U^ on M) by the flows (2.17) and (2.18). 
(Here /c" are m quasiperiods of the function ^{x).) So we obtain that our 
family of m-phase solutions of (1.6) is also the family of m-phase solutions 
for (1.7) and (2.11) and we assume also the existence of the solutions (2.2) 
for these systems based on the family M.. 

We can also conclude that in our situation the variational derivatives of 
the functionals (2.15) and (2.16) with respect to ^{9) at the points of the 
submanifold M. are the linear combinations of the corresponding "right eigen 
vectors" K(^q){9 + 9q,U) (see [7]- [10] and references therein). Indeed from 
the conservation of the functionals (2.15) and (2.16) by the flows ip\a = ^\a 
and 

(^j = g^(v9,A;Ve-,...) 

we can conclude that the convolution of their derivatives with respect to v^*(^) 
with the system of constraints (2.9) is identically zero for all the periodic 
functions with respect to all 6'" and for any k^, . . . ,k^ and uj^, . . . , lo^ . So 
we can take the variational derivative of the corresponding expression with 
respect to <f-'{9') and then omit the second variational derivative of I'^ and H 
according to the conditions (2.9). After that we obtain that the variational 
derivatives of J*^ and H are also orthogonal to the image of the linearized 
operator (2.9) at the points of Ai and so are the linear combinations of 
i^{g){0 + 0o,U) onM. 

Lemma 2.1 

Suppose we have the properties (A)-(C) and (I)- (II) for our family of 
m-phase solutions of (1.6). We put 

U" = (V) = F (2.19) 

on the space M. and then define the functions k" = k°'{U) on the submanifold 
M. 

Then the functionals K°' = k°'{I[ip]) on the space of 271 -periodic with re- 
spect to each 9°" functions (and also at the space of quasiperiodic functions 



ip{x) with m quasiperiods) have the zero variational derivatives on the sub- 
manifold A4 . 

Proof. 

As we have from (II) the maximal number of the linearly independent 
variational derivatives oi I^ on M. is h + m = N — m. So we have m linearly 
independent relations 

N .fu 

T.KiU)—-^0, a = l,...,m (2.20) 

{(p{9) = {(p^{9), . . . , ip"^{9))) being considered at the given k^, . . . ,k"^ at any 
point of Ai (or in other words 

N ^ju 

^A^(f/)--- = 0, a = l,...,m (2.21) 

when considered in the space of functions with m quasiperiods.) We have 
here the standard expression for the variational derivative and the formula 
(2.14) for I". 

Since we can obtain the change of these linear combinations of the func- 
tionals I" on A1 only due to the variations of k in (2.14) but not of <f^{0) (or 
in other words only if we have the non-bounded variations of (f^{x) after the 
variations of the quasiperiods) we have on Ai 

N m 

Y.K{U)dU'^ = Y.4"\U)dk^iU) (2.22) 

for some functions fi^ (U). 

If W^ are the coordinates on A4 then the matrix /il" has the full rank and 
is reversible. So we obtain the differentials dk^ as the linear combinations 
of differentials J2u=i ^^{U)dU^ corresponding to the functionals with zero 
derivatives on M. 

m N 

dkf' = Y.if^-%^Y.^^^\u)du^ 

So Lemma 2.1 now follows from (2.20). 
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Remark 1. 

As can be seen from the proof of Lemma 2.1 the variational derivatives of 
I^ on M. should span all the (A^— m)-dimensional linear space generated by all 
f^[q){0 + 9Q^ U) if we want to take {V^) as the independent coordinates on M.. 
It is essential also that we consider the full family of m-phase solutions given 
by (2.9) at different u and k but not its "subspace" and have m independent 
relations (2.22) on N differentials dU^ from m relations (2.20). 

Remark 2. 

We note here that the equations (2.21) were introduced at first by S.P. 
Novikov in [11] as the definition of the m-phase solutions for the KdV equa- 
tion. 

We now prove a technical lemma which we shall need later. 

Lemma 2.2 

We introduce the additional densities 



nr(.)(y,y.,...)^ nT""^ (2-23) 



for k > 0, where ip\,j. = d^ / dx'^ ip'' {x) . 

Then on the submanifold M. we have the relation 



^ f)k°' 1 /"Stt /•27r 

,Ss7^(2^/„ I i:P*'^-(£/).^^*'^'-(£/)*;.„,«.,^,„..(».t/)x 

xnr(ri(<l'(fc.)(»,C),*^<l'(„.)<)-,(»,C),...)rf"'« = S"^ (2.24) 

at any U and 9q. 

Proof. 

According to Lemma 2.1 we should not take into account the variations 
of the form of $(«„) {0 + 6q, U) when we consider the change of the functionals 
A;"(/) on M.. So the only source for the change of these functionals on M. is 
the dependence on the wave numbers k in the expressions 
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27r /-Stt 



(27r)™7o '"Jo 
So we can write 



V•'(i^„.„k-•t^^.„fe-„ Ffc*<l.,j„,„,„,, . . .)rf"'9 



where the values of dl^[ip\/dk^ on M. are given by the integral expressions 
from (2.24). Since the values of the functionals A;°(/) on M. coincide by the 
definition with the wave numbers /c" we obtain the relation (2.24). 

Lemma 2.2 is proved. 

For the evolution of the densities V'^{ip, ^Px, ■ ■ •) according to our system 
(1.6) we can also write the relations 

jP-'i^, ^x, . . .) = OxQ-'^'i^, ^., . . .) (2.25) 

and Whitham's system (2.3) can be also written in the following "conserva- 
tive" form 

OtU'' = dxiQ"'') z/ = l,...,iV (2.26) 

for the parameters W^. 

From the existence of the series (2.2) on the space of 27r-periodic with 
respect to 6 functions it can be shown that this form gives us the system 
equivalent to (2.3) in the generic situation and all the other local conservation 
laws of form (2.25) (if they exist) give us the equations 

dT{V) = dx{Q'') 

compatible with (2.26) for the full set of parameters W^ (see [1]- [10]). 

The conservative form (2.26) of Whitham's system will be very convenient 
in our considerations on the averaging of Hamiltonian structures. 

Now we make some "regularity" requirements about the joint Ai of the 
submanifolds A4uj,k at all u and k corresponding to the m-phase quasiperiodic 
solutions of the system (1.6) i.e. 
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(Ill) We require that M. is an {N + m)- dimensional submanifold in the 
space of the 27i-periodic with respect to each 6" functions. 

The property (III) means nothing but the fact that the shapes of the 
solutions of (2.9) are all different at different uo and k in the space of the 
27r-periodic vector-functions of 9 so that oo and k can be reconstructed from 
them. It is easy to see that this requirement corresponds to the generic 
situation. We shall use the property (III) for our procedure of averaging of 
bracket (1.1). 

We shall work with the full family of 27r-periodic solutions of (2.9) the 
functions of which will also depend on the "slow" variables X and T. To 
define this submanifold in the space of functions (/9(^, X, T) with the system of 
constraints like (2.9) we should extend the coordinates V^ as the functionals 
of ^{0) in the vicinity of our submanifold M.. This can be easily done (see 
[14]) by the following way: 

Let introduce A^ different functionals 

such that their values A'^ are functionally independent on the functions from 
the submanifold M.. Then we can express U^ = f/^(A) in terms of A" on 
A^ and after that extend them as the functionals [/"{A) on the space of 
27r-periodic with respect to each O"' functions. 

We shall also expand the coordinates 6q (see [14]) in the vicinity of Ai 
by introduction of, say, functionals 

which are zero if ^{9) = <I>(j„)(^, [/[(/?]). In the generic situation we can 
locally express the values of 6q on Ai in terms of 5" and after that put 
9q = ^Q (5 [(/?]) in the corresponding local coordinate maps in the vicinity of 
M. 

Now we consider the system 

V\e, X) - $J,„)(^ + OMMV]) = (2.27) 
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where 6q[^p] and If^lip] are the functionals in the vicinity oi Ai as the system 
of constraints which defines Ai in the space of 27r-periodic with respect to 
each 6°" functions. 

We can see now that the hnearized system (2.27) 

I /■2-K /■2-K 



{2tt)'^ Jo "Jo 
where 



L][u,eo]((^,e')6<^^{e')]d-e' = o 



L)yuM^e,e')^5]5{9-9')- 



5e^\^] A.. 5f/o"[ 



Y,^l^^^^^{e+e,[^]Mv])^Y^^-T.^\^nM^+Hv]Mv])^ 



a=l 



SV'{0') ^1 ^"^ UL7^J,^L7^J/ ^^j^Q,^ 



at any point {U,9q) of A4 has exactly N + m solutions C,{p)[u,eo]{^) corre- 
sponding to the tangential to M. vectors ^(in)e" and ^(in)^, a = 1, . . . ,m, 
u=l,...,N. 

It is evident also that all the "left eigen vectors" K{p)[u,0o]{(^) orthogonal 
to the image of L are given by the variational derivatives 66Q[(f>]/6if>'{6) and 

From the invariance of the submanifold Ai with respect to the fiows (1.7) 
and (2.11) we can also write here the relations 

I f2w r2n 

-^}[t/A)](^'^')x 



(27r)'" Jo "Jo 

X 5^,) ($(,„) (^' + 9o, U), A;"<l>(,„)e.(^' + 9o, U), . . .)dre' = (2.28) 

and 

2w r2w 



{2ti)^Jo "'Jo "^■[^•^"1 

xgj,)($M(^' + ^0, f/), A;"<l'(.„)e4^' + ^0, f/), . . W'O' = (2.29) 

for any i,k at any point {U,6o) of 7W where k" = A;"[$] can be considered 
now as the values of the corresponding functionals on A^ . 
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We now introduce the space of functions ip{9,X,T) depending on "slow" 
parameters X and T and 27r-periodic with respect to each 6°'. The systems 
(2.27) being considered independently at different X give us the system of 
constraints defining the subspace Ai' in the space of functions ip{6,X) cor- 
responding to m-phase solutions of (1.6) depending also on the parameters 
X and T. 

Now after the introduction of the "modified" constraints (2.27) 



Glu,e„] {0, X) = —— ... L^yjj^s^^ {6, 6') x 



1 /'27r /'27r 

(27r)™ io "Jo 

X {^^{6') - %^^{e' + ^oM, f/M)) d-^e' (2.30) 

we can take W^^X), 9q{X) and G^rjjgi{9,X) such that 
^f...f .,,.j(. + .o(X))x 

xq^e,^](^,X)rf'"^ = , p=l,...,N + m (2.31) 

as the coordinates in the vicinity of Ai' instead of the <^*(^, X) since we can 
find uniquely ip^{6,X) from 

(2^ /" • r^PM (^' ^') {^' (^') - %n) iO' + 0o,U)) d-9' = q^,,,, (9, X) 



if we have the conditions (2.31) and the values of f/'^(X) and 9q{X). 
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Remark. 

Certainly we have here some freedom in the choice of the constraints 
G^{9,X). For example we can take also the expressions (2.27) as the system 
of constraints defining Ai' . We prefer here to take the constraints in the form 
(2.30) just to fix the uniform orthogonality conditions (2.31) in the vicinity 
ofA^'. 



We shall need also another coordinate system in the vicinity of Ai' which 
differs from the described above by the transformation depending on the 

^This system of constraints is different from the system introduced in [14]. 
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small parameter e and singular at e — » 0. Namely, we recall our integrals 
(2.5) 

F = Jv''iv,V,,...)dx 
make a transformation X = ex and define the functionals 



r{x) = ^ 1^'^. . . j^^ v^ivie, X), e^xie, x), . . .)rf"^ (2.32) 

on the space of 27r-periodic with respect to each 6°' functions ip{6,X). 
Let us also introduce the functionals 



C(^) = ^0 (^) - ^0 (^o) -- r r ( J(X'))rfX' (2.33) 



1 f^ 
for some fixed point Xq. We have identically 



^*"(Xo)=0 (2.34) 

As was shown in [14] we can also obtain the values of If^lX) and Oq{X) 
from J^{X), ^q"(X) and 6q{Xq) on Ai' as the formal series in powers of e 
and we shall have for these series 

W'ixu e*,] = r{x) + Yl ^'u^.p, Jx, e^x, ■ ■ •) (2.35) 

k>l 



e^{x)[j,e*] = e*^{x) + 0^{Xo) + - ['' k'^iJix'))dx' (2.36) 

e JXo 

The form of the relation (2.35) will be important in our considerations, 
so we reproduce here the calculations from [14]. 

We remind that the values J'^iX), O^'^iX), e^{Xo) and f/^(X) are con- 
nected on Ai' by the relations (the definition of J'^(X)): 

r{X) = ^ p . . fj V^i^Ud + s{X^ e), f/(X)), 

e^x<^i^n){e + s(X, e), U{X)), . . .)rf™^ = 
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'''V^($(,„)(^ + s(X,e),f/(X)), 



(27r)™ Jo "Jo 

A;"( J)9e.$(.„)(^ + s{X, e), U{X)), . . .)(r'e+ 

+ ^^e'^[^...jyU^Ue + s{x^e)Mx))....)d-e 

where 

s(X, e) = ei{X) + ^o(Xo) + - /'' k{J{X'))dX' 

and Pn;)($(m)(6' + s(X, e), . . .) are local densities depending upon $(in)(6' + 
s{X, e), U{X)) and their derivatives with respect to U" and 6" with the coeffi- 
cients of type: UxiX),UxxiX), . . ., k{J),dxk{J),d\k{J), . . ., and 
6lx{.X),6lxxi.^)^ ■ ■ -1 given by the collecting together these terms, having 
the general multiplier e^ . The term corresponding to the zero power of e is 
written separately. 

After the integration with respect to ^, which removes the singular at 
e — !► phase shift 6q in the argument of $(«„), we obtain on Jvl': 

J^iX) = ("{J, U) + ^e C(k){U, Ux, . . . , Ukx, J, Jx, ■ ■ ■ , Jkx, ^ox; • • • ; ^okx) 

k>l 

(2.37) 
The sum in (2.37) contains the finite number of terms. The functions Co;) 
and C^ are the integrated with respect to 6 functions 'P^^^ and V^ respectively. 
So, since 



CUU) = ^|J:.^V(<i>(.„)(^,f/),r(j)<i>(,„),.(^,f/),...)rf'"^ 

we obtain that the system 

r{X) = C{J{X),U{X)) (2.38) 

is satisfied by the solution J'^{X) = [/"{X) according to the definition of 
the parameters If^. Since we suppose that the system (2.38) is of the generic 
form we shall assume that (locally) this is the only solution and put J'^{X) = 
If^lX) in the zero order of e. 
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After that we can resolve the system (2.37) by the iterations taking on 
the initial step U^iX) = .r{X). The substitution of (2.35) into (2.37) under 
the condition of the non-singularity of matrix || qjjI \\\u=j will sequentially 
define the functions -Un,-). So we obtain the relations (2.35) and (2.36). 

Now we can take also the values of J^{X)^ 91°" [X]^ 9q{Xq) and 
G?^r^i g^r^ii(^,X) with the restrictions (2.34) and also 

(2^ C- ■ ■ r '^(^)[^[^i(^)] (^ + ^o(^) + ^°^^°) + 7 ixo ^('^(^'))^^') ^ 

as the coordinates in the vicinity of Ai'. 

We define now the Poisson bracket on the space of functions ip{6,X) by 
the formula 

Wie, X), ^^{9', Y)} = Y: Bl\^, e^x, . . .)e'5^'\X - Y)5{d - 9') + 

fc>0 

+ -5{9 - 9') Y: e,Sl,)i^, eifx, ■ ■ -MX - V')5J,)(^, eify, . . .) (2.40) 

^ fe>0 

which is just the rescaling of the (1.4) multiplied by 6{9 — 9'). We normalize 
here the 5-function 5{9 - 9') by (27r)". 

The pairwise Poisson brackets of the constraints G\^g A9^X) on M.' can 
be written in the form 



{G\9,XIG\9\Y)]\m> = (^/" ••/"^M^™o(x)](^,r)x 

x^i[i/(y)A,(r)](^» X W\r,X),^%a,Y)}\M'd'-rd-a (2.41) 

(we can omit the Poisson brackets of L\ and LionAi' and also the brackets of 
the functionals 6*0 [y?] and W^l^p] from <l>(j„) in (2.30) since they are multiplied 
by the convolutions of the corresponding L-operators with the "right eigen 
vectors" $(m)e« and ^(in)U'^ which are zero on Ai'). 
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Brackets (2.41) evidently satisfy the orthogonality conditions 



wr r • • • r"'^('^)^[^[^'^oi(^)] (^ + ^o(^) + ^o(xo) + \ £ k{j{x'))dx^x 

x{G\e,x),G^{e',Y)}\M'<r^d = o (2.42) 

''\G\e,x),G^ie',Y)}\M'y< 



(27r)™yo '"Jo 

xk^.muijmy)] (^ + ^o(Y) + ^o{Xo) + -J'^ k{j{Y'))dY'j <rd' = (2.43) 

for q = 1, . . . ,N + m in the coordinates J{X), 0q{X) and ^o(-^o) on the 
submanifold Ai'. 

We note now that every derivative with respect to X or F appears in the 
bracket (2.40) with the multiplier e but being applied to the functions 



^'\e, X) = $J,„) (^e + e*,{x) + ^o(Xo) + ^ /^"^ k{j{x'))dx', u[j, e;]{x)^ 

(2.44) 
on Ai' contains the nonzero at e ^ term k°'{J)d/d6°'. Now we formu- 
late the statement about the structure of the bracket (2.40) on Ai' in the 
coordinates J{X), 6q{X) and 6o{Xq). 

Lemma 2.3 

The pairwise Poisson brackets of constraints G^jjq ^ {6, X) have no singular 
terms at e —^ and have no non-local terms in the zero order of e (e*^) on 
M' at any fixed coordinates J^iX), 6'*"(X) and 9'^{Xq) (such that U{X) = 

u[j,e*]ix), doix) = eo[j,e*„eoiXo)]ix)). 

Proof. 

The first statement is evident for the local part of bracket (2.40) since 
any differentiation with respect to X in it appears with the multiplier e and 
has the regular at e ^ form k°'{J{X))d/d9°' + 0(e) being applied to the 
functions of form (2.44). So we should check only the non-local part of (2.40) 
which contains the multiplier e~^ in it. But according to the relation (2.28) 
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and also (2.35) we have that the terms arising on the both sides of z/(X — Y) 
(the convolutions of L with S(^k){^, k^'^ga, . . .)) are of order of e on the Ai' 
in the coordinates J^iX) and 6q°'{X). So we obtain that all the non-local 
part of (2.41) is of order of e on Ai' at any fixed coordinates J'^(X), ^q"(X) 
and ^^(Xo). 

Lemma 2.3 is proved. 

We formulate now the last "regularity" property of the submanifold AX' 
with respect to the Poisson structure (2.40). 

We consider in the coordinates J'^(X), 91°'{X) and ^o(-^o) on M' the 
linear non-homogeneous system on the functions 

f^m] (o' + eiiY) + ^o(Xo) + i f^^ k{J{Y'))dY', Y, e) 
having the form 

T /'27r /'27r 

J2^ Jo "Jo ^'~^[^M'^oM](^'^)''^[t/M,eoM](^''^)}l-^'^ 
xfj U + ^ol^') + MXo) + \f^ k{J{Y'))dY\ F, e J d^'O'dY = 

= {Glu^^]^SoM]ie,X),F[^]{e)}\M' (2.45) 

where F[(/9](e) is the functional defined in the vicinity of Ai'. 

After all differentiations with respect to X we can omit the term ^o(^) + 
Oo{Xo) + ^ /x,j k°'{J{X'))dX' which occurs in all functions depending on 9 
and X in (2.45) and then consider the system (2.45) at the zero order of e. 

From Lemma 2.3 we have that at the zero order of e the brackets 
{G'^{6,X),G^{6',Y)} on At' do not contain the non-local terms depending 
on z/(X — Y). For the derivatives with respect to X which arise with the 
multiplier e from the local terms of {ip''{T,X),ip^{(r,Y)}\_M' we should take 
in the zero order of e only the main part k°^{J)d/d9'^ . So in the zero order 
of e we obtain from (2.45) just the linear systems of the integral-differential 
equations with respect to 9 and 9' on the functions fj{9',X) independent at 
different X. We have also that the right-hand side of (2.45) satisfies at any 
X and e the compatibility conditions (2.42) (let us remind that U'^[J,9q] are 
the asymptotic series at e ^ 0). 
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(IV) So we require that the system (2.45) is resolvable on Ai' for any 
F[ip]{e) in the class of 27r -periodic with respect to all 9°" functions and its 
solutions can he represented in the form of regular at e ^ asymptotic series 

kAk) 



fn.mi9,Y,e) = j:e''f^ll,.^ie,Y) 



n>0 



for regular at e -^ right-hand sides of (2.45). 

The condition (IV) is responsible for the Dirac restriction of the bracket 
(2.40) on the submanifold A^'. 

Now we prove the statement which will be very important for our aver- 
aging procedure. 

Lemma 2.4 

The Poisson brackets of the functionals ^q"(X) with J'^iY) are of order 
of e at e ^ on the M! at any fixed coordinates J^[X) 6'q"(X) and Q'^{X^ 

{ei\X), r{J)}\M' = 0(e), e ^ (2.46) 



Proof. 

First we note that the Poisson brackets of ip'^{9,X) with the functionals 
J'^{Y) can be written in the form 



W{e, X), .r{Y)} = Y: Cn^ie, X), e^O, X),.. .)e'5^'\X - Y) + 

fc>0 

+ J2^'^Sl,^i^i(^,X),e^x{e,X),...)u{X-Y){Fl',^i^{e,Y),eM0,Y),...: 



'Y 
fc>0 

for some C]^{'{>, ^^x-, ■ ■ •) and F^^^i^p, e(/9y, . . .) (we have integrated with re- 
spect to ^'). 

So the flow generated by the functional / q{Y^J^{Y^dY (where q{Y^ has 
a compact support) can be written as 

fc>0 
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+ E efc^w(v^, e^x, ...) jv{X- Y)q{Y)^F(,^{^, e^y, . . .)dY = 
k>0 '' "'^ 

= J2Cl^{v,eipx,..yqkx{X) + 

fc>0 
fe>0 

- E e.^(\)(v^, eifx, ...)jy{X- Y)F^,^i^, eify, . . .)qY{Y)dY (2.47) 



A;>0 



As can be easily seen the local terms of (2.47) have the form 



q{x) 



Cq{v, ev^x, • • •) + E efe'5(fe)(v^, ev^x, • • •)^(fc)(v^> ey^x, • • • 



fe>0 



+ 0{e) 



where the term in the brackets is just the flow generated by the functional 



1 f fiiv fiiv 



(27r 

In the non-local part of (2.47) (the last expression) we have the convolu- 
tion of the "slow" functions qviX) with the rapidly oscillating FY^Aip, ey^y, . . .) 
where ip'^{d^Y) has the form (2.44). So in the leading order of e we can ne- 
glect the dependence on 9 of the last integral of (2.47) and take the averaged 
with respect to 6 values {Fy^\) on Jvi' instead of the exact FK^Jip, eipy, . . .) in 
the integral expression in (2.47). 

After that we obtain that the non-local term of (2.47) gives us in the zero 
order of e the linear combination of the flows 5'(fe)((/9, eipx, ■ ■ •) considered on 
the functions 



ip\e, X) = %^) (e + e*{x) + ^o(Xo) + - f" k{j{x'))dx', j{x) 



e Jxo 



at any fixed point X. 

From the invariance of the submanifold A^ with respect to the flows (2.17) 
and (2.18) we can conclude now that the flow (2.47) being considered at the 
points of Ai' with flxed coordinates J{X), 6q{X), 6q{Xo) in the zero order 
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of e leaves Ai' invariant and generates on it the linear evolution of the initial 
phases 

e^ix) = e*-ix) + e^ix,) + 1 T k^iJix'))dx' 

with some frequencies Q?^{X) . Here we use the formula (2.35) for U[J,6q] 
and we can claim now that the Poisson brackets of the functionals Oq{X) with 
Jq{Y)r{Y)dY at the points of M' with fixed coordinates J^iX), Ol'^^X) 
and ^o(-^o) have the form 

{e^,{X), j q{Y)r{Y)dY} = n'(:^[J,ei]{X) + 0{e) (2.48) 

We now prove the relation 



{k^{J{X)), j qiY)riY)dY} = e ^n'(:^[iei]{X) + 0{^) (2.49) 

at the points of M' with the fixed values of J''(X), ^o"(^) and ^o (Xq). 

Using again the relation (2.35) we can write for (2.47) at the points of 
M' 



^\ = nf;^{x)^i,^^,,(9 + e*{x) + ^o(Xo) + -J\{j{x'))dx', u[j, e*]{x)j+ 
+ev' (e + e*{x) + ^o(Xo) + i J^ k{j{x'))dx', [j, e*,]\ (2.50) 

where [J, 9q] means the regular at e -^ dependence on J, Jx, Oqx, ■ ■ ■■ 
We are interested in the evolution of the functionals 

We have 
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r-27r /■27r 



//•••// (nf(o)V^i + ^tiiMx + nf(2)6Vixx + ■ ■ ■) d-e 



(2vr) 

where the densities H^.^n were introduced in (2.23). 

It is easy to see that (2.50) does not change J^{X) at the zero order of e 
and we can also state that the terms of order of e in (2.50) (i.e. eri'^{6+. . . , X)) 
are unessential for the evolution of k{J{X)) on Ai' at the order of e. Indeed, 
their contribution to the evolution of J'^(X) in the order of e is 



^ (^ C ■ ■ r K)V' + ll':afVx + ^^../rf^x + . . .) d-e (2.51) 

where we should take only the main term k"'{J{X))d/d9"' for the derivatives 
ed/dX in the formula (2.51). After the integration by parts we have for this 
contribution 



'Wr C- ■ ■ C ("•'»' ' *''^"''- ^ ■ ■ ■) "'<" ^ ■ ■ ■ • ''""''' 



But after the substitution of the main part of v^*(^, X) 

^U [o + ^o(^) + ^o(Xo) + \ J^ k{J{X'))dX', J{X) 



(according to (2.35)) into the densities Il^i^Jip,eipx, ■ ■ ■) we obtain in the 
leading order of e the convolution oiri{9, X) with the variational derivative of 
the functional /^ introduced in (2.15) with respect to ip{6, X). Our statement 
follows now from Lemma 2.1 which claims that the variational derivatives of 
the functionals A;" (J [(/?]) are identically zero on the space of m-phase solutions 
of (1.6). 

Consider now the first term of (2.50). We have that the evolution of J^{X) 
which is responsible for the evolution oi k{J) is given by the expression 

^J-(X)=fif^^(X)x 

1 j'2-K j'2-K JjT>l^ 
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x<^l,^^,,{e + s{X,e),U[J,e*]{X))d^-e + ©(e^) 

where s(X, e) = 9*{X) + 9o{Xo) + i j% k{J{X'))dX' . 

The first term here after the substitution of exact (p* in the form 

if\e,X) = ^l^^^,,i9 + siX,e),U[J,9*,]iX)) 
on the Ai' as can be easily seen is just 

while the second term on Ai' in the leading order of e is equal to 






^ X 
p>\ 

xnf(^) ^$(,„)(^ + s(X, 6), J(X)), F( J(X)) A$(,„)(^ + s(X, 6), J(X)), . . .jx 

xrHJ(X)) . . . A;"-H^(^))<^U/^.^i...r.-i (^ + s(X, e), J(X)), . . .) rf™^ 

(2.52) 
which coincides with the integral expression from (2.24) in Lemma 2.2. So 
from Lemma 2.2 we have that the summation of (2.52) with dk'^/dJ^ is equal 
to e (n^^JX)] 5^ and we obtain that 

^tnJ)=e^nf^^{X) + 0{e') 

i.e. the relation (2.49). 

Now using (2.48) and (2.49) we can write that 

{9;^iX),JqiY)riY)dY} = 

1 r^ 

'xo 



{9^{X) - 9^{Xo) -- r nJ{X'))dX', fq{Y)rdY} = 0{e) 
e Jxq J 
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for any q{Y) in our coordinates on Ai' . 
So we have 

{e*,-{x),riY)}\M' = o{e) 

at any fixed coordinates J^(X), O^'^iX) and 9^{Xo) 



Lemma 2.4 is proved. 



3 Averaging procedure. 

Now we formulate the averaging procedure of the Poisson bracket (1.4) on 
the family of m-phase solutions of (1.6) under the conditions of "regularity" 
formulated above. 

Theorem 3.1 

Suppose we have the Poisson brackets (1.4) (ind the local system (1.6) 
generated by the local Hamiltonian function 



H = I VHiv,<^x,---)dx 



which has N{> 2m) -parametric full family of m-phase solutions modulo m 
initial phase shifts 9q . 

Let us have also N commutative local translation invariant integrals 



V''{(p,(p^,...)dx 

{F,H} = 0, {r,F} = o 

which generate local flows according to the Poisson bracket (1.4) o-nd the 
averaged densities of which can be taken as a parameters U" on the space of 
m-phase solutions of (1.6) (the conditions (A)-(C)). 

Then under the "regularity" conditions (I)-(IV) for the space of m-phase 
solutions of (1.6) we can construct a Poisson bracket of Ferapontov type (1.5) 
for the "slow" parameters f/'^(X) which coincide with the averaged densities 
{V''){X) by the following procedure: 

We calculate the pairwise Poisson brackets ofV'^{ip, fx, ■ ■ ■) in the form 
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{V'iv, ^., . . .),v^{v, V.,, ...)} = E ^7iv, v., ■ ■ ■)S^'\x - y)+ 



fc>0 



+ E efc te)(<^, ^x, ■ ■ -n^J^ix - y) (FL{ip, ify, . . .) 



fc>0 



(where are finite numbers of terms in the both sums). Here we have the 
total derivatives of the functions F^j^-. and Fn,\ with respect to x and y as a 
corollary of the fact that both I" and I^ generate the local flows (2.11). From 
the commutativity of the set {f^} we have also 



(3.1) 



fc>0 



for some functions Q^^l^p, ipx, • • •)■ 

Then for the "slow" coordinates If^lX) = {'P^){X) we can define the 
Poisson bracket by the formula 



{f/'^(x),f/^(r)} 



{AT){x) + E e. mu)Fi:,;) - (i^r.))^^^) w 



fe>0 



5'{X - Y)+ 



9(g'^^)(x) 



dX 



J2^k 
k>0 



5(^r.))w 



dX 



(^S))W 



6{X-Y)- 



j:,,!<£|P.,A--r).^<^«)<^> 



fc>0 



dX 



dY 



(3.2) 



where the averaged values are the functions of U{X) and U{Y) at the corre- 
sponding X and Y . 

Bracket (3.2) satisfies the Jacobi identity and is invariant with respect to 
the choice of the set {F, . . . , l'^} satisfying (A)-(C) if it is possible to choose 
such integrals in different ways, i.e. 

ifW = (V), U^ = {V") and {U''{X),U''{Y)}, {U''{X),U>'{Y)}' are the 
brackets (3.2) constructed with the sets {I'^} and {/^} respectively then 
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{U'^{X),U^{Y)y^^X){U\X),U^iY)}^Y) 

Proof. 

The most difficult part is to prove the Jacobi identity for the bracket 
(3.2). For this we use the Dirac restriction of the Poisson bracket (2.40) 
on the submanifold Ai' with the coordinates J'^(X), ^q"(X) and 6q{Xo) on 
it. According to the Dirac restriction procedure we should find for J''[X), 
e*Q''{X) and e^{Xo) the additions of type 



vix) = j^J'J..fJ jv';[iei'^,eo{x,)]{x,6\z)&{6\z)<r^e'dz 



and 



o" = ^ p. . fj J o'^iJ, c, 0o{XoW', z)G^{e', z)d^e'dz 

such that the fluxes generated view (2.40) by the "functionals" J'^{X) + 
ViX), ^(^"(X) + W{X) and e^{Xo) + O" leave M' invariant i.e. 



^ ^'" ''^\g\9, Y), G^{9', Z)}\m' X v^'JX, 9', Z)d"'9'dZ 



(27:)'^ Jo "Jo 

= -{G\9,Y),r{X)}\M' (3.3) 



-^ j\ . . fj{G\9, Y), G^{9', Z)}\m' X w^{X, 9', Z)d"'9'dZ 



^{G\9,Y),9*'^{X)}\m' (3.4) 

37 



^ ^'" ''^\G\e, r), G\e', z)}\m' x o"(^', z)(re'dz 



(27r)™yo ■■'Vo 

= -{G^(^,F),^o"(^o)}U' (3.5) 

and put after that for the Dirac restriction on Ai' 

{ j^(x), r{Y)}D = {r{x) + vix), r{Y) + v^{y)}\m' = 
= {r{x),j^{Y)}\M' - j^ fj. ■ ■ fj v'({x,e,z) xv^(Y,e\z')x 
x{G\e, z), G^{e', z')}\M'cre(re'dzdz' (3.6) 

and by the same way 

{.r{x),er{Y)}o = {r{x),e*^iY)}\M'- 

I r2TT I-2TT 



v'^{x,e,z)xwUY,e\z')x 



(27r)2™yo ••Vo *' ' ' ' ^ 

x{G\e,z),G^{e',z')}\M'd"'ed"'e'dZdz' (3.7) 

{enx),ef{Y)}r, = {er{x),e;f'{Y)}\M'- 

I r2n r2TT 



w"{X,e,Z) xwUY,9',Z')x 



(27r)2™yo "Jo ^ 

x{G\e, z), G^{e', z')}\M'd"'ed"^e'dZdz' (3.8) 

and so on. 

After the calculation of the brackets in (3.3)-(3.5) and the substitution 
of ip{9, X) in the form (2.44) we obtain the regular at e ^ systems for the 
functions 'yj(-^, 0, Z, e), w'j{X, 9, Z, e) and o'j{9, Z, e) such that 

v'^iX, e, Z, e) = t;J (x, 6 + 0*(Z) + ^0(^0) + ^ ^^ k{J{Z'))dZ', ^ 
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wJ(X, e, Z, e) = wj (x, 9 + 91{Z) + ^0(^0) + - /^ KJ{Z'))dZ' 

\ t JXn 



and 



one, Z,e) = o'^{e + e*{Z) + ^o(Xo) + - [^ k{J{Z'))dZ', e 



Xa 



which coincide with the system (2.45). 

From the arguments analogous to those used in Lemma 2.3 and the fact 
that the flows generated by the functionals J^{X) leave invariant the sub- 
manifold A^' at the zero order of e in the coordinates J{X), 61{X), ^o(-^o) 
we have also that the right-hand sides of these systems are regular at e — >■ 
in these coordinates. 

So according to (IV) we can flnd the functions -uj, w" and o" in the 
form of regular at e ^ asymptotic series, (the functions fJ(X, 6*, Z, e), 
w"{X,6,Z,e) and o'j{9,Z,e) are not uniquely deflned but it can be shown 
that it does not influence on the Dirac restriction of the bracket (2.40) on 
Ai' according to the formulas (3.6)-(3.8)). 

Besides that, as was mentioned above, the flows (2.47) generated by the 
functionals J q{X)J^{X)dX on the functions of form (2.44) leave invariant 
the submanifold Ai' at the zero order of e and generate the linear evolution 
of the initial phases. So we can conclude that the right-hand side of the linear 
system (3.3) contains no zero powers of e and we should start the expansion 
for v^X, 9, Z, e) from the flrst power. 

Now we have 



v-^lJ, 91 ^o(Xo)](X, 0. Z,e) = Y. e'^J(fe)[J, ^o> ^o(Xo)](X, ^, Z) 

k>l 



w'^lJ, 91 ^o(Xo)](X, 9, Z,e) = Y. e^w^^^iJ, 0^ ^o(Xo)](X, 9, Z) 

k>0 

o'^[J,9*„9o{Xo)]{9,Z,e) = Y.e'o'^^,)[J,9*,9o{Xo)]{9,Z) 

k>0 
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According to the relations above and Lemma 2.3 we can see now that the 
corrections to the values WiX), J^(F)}|^/ and {^o"(^)> J^iY)}\M' due to 
the Dirac restriction on Ai' are of order of O(e^) and 0(e) respectively. 

According to the relation (2.35) we can also substitute the values J'^{X) 
instead of f/''[J, ^q](X) in the leading order of e as the arguments of the 
averaged functions on Ai' 

Now we calculate the values {J^(X), J^^(Y)}\_m'. 

We consider the Poisson brackets 

{j q'{X)r{X)dX, j q^{Y)J^(Y)dY}\M' (3.9) 

for the arbitrary smooth q'^{X) and q^^(Y) with compact supports. For the 
Poisson brackets of the densities V^{ip, eipx, ■ ■ •) according to (2.40) we have 
the expression: 



{v^{e, x),v'^i9', y)} = y: ^r(^, e^x, . . y6^'\x - Y)6{e - e')+ 

k>0 



+e6{e - e') Y: eu (F^,){v, e^x, ■ ■ .)) ^ v{X - Y) (F^^^i^, e^y, . ,^^ 

A;>0 

such that 

{r{x),j^{Y)} = 
= E ^'t:^ r ■ r Ai^{v{e,x),evx{e,x), . . .)<re 6^'\x - y)+ 

yv(X-Y) X (F^,(f,(e,r),epv(O,y),^^0)^ (310) 

Now we should substitute the functions v^*(6', X), ip\6,Y) on Jvi' in the 
form 



$5,„) [e + ei{x) + ^o(Xo) + \ j^ k{j{x'))dx', u[j, e*]{x)\ (3.11) 
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and 



*M {0 + 0*oiy) + UXo) + \ £ k{J{r))dY', U[J, e*]{Y)^ (3.12) 

respectively. 

It is easy to see that the local part of (3.10) gives us the expression 

1 f2TT f2TT 



(27r)™ Jo "Jo 

_d_ 
dX 



xA^,''(<i>l,^^{e + s{x),u[j,e*]{x)),e^<^l,^){e + s{x),u[j,e*]{x)),..] 



x6{X-Y) + e{A''/){J{X))6'{X-Y) + O(e^) (3.13) 

in the coordinates J(X), 9^{X) and 6^0 (Xq) on A^' where s{X) = 9^{X) + 
eo{X,) + \!lk{J{X'))dX'. 

Here we used only the main part of (3.11) $(j„)(6' + s{X), J{X)) in the 
second term of (3.13) and replaced U"{X) by J^i^X) according to (2.35) in 
the arguments of the averaged functions modulo the higher orders of e. 

In the non-local part of (3.10) we have for (3.9) the following equalities: 



1 f2. ^2.1^ ,^^^^^gF(-)(«i>(,„)(g + g(x), [/(x)),...; 

X 



(27r)™Jo Jo e^o 



{2^rJo--Jo e^,^^- -- dX 



xz/(X - F) X eg^(F) — ^-^^ — d 6 



J J (27r)™ io io ^0 (9X(9F 

xF('^,) (<i.(,„)(^ + s{x), u{x)), . . .) Fii^ (<i.(,„)(^ + 5(r), t/(r)), . . .) d'-e 



dXdY-1- r. .. Tj: e, 6[g^(X)z/(X - Y)q^iY) + 

[ZTV) JO JO j^>Q 
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+ (g'^(X)g^(F) - q\{X)q^{Y))6{X -Y)- q^{X)q'^{Y)6'{X - Y)] x 

xFi',^ ($(,„)(^ + s{x), u{x)), . . .) Fii^ (<i.(,„)(^ + s{Y), u{Y)), . . .) <re = 

xF('^,) ($(,„)(^ + s(X), U{X)), . . .) FJ) (<l>(,„)(^ + siY), UiY)), . . .) d"^9+ 
+e E e'^ / iQ''iX)q^xiX) - q^xiX)q^iX)) {Fi:,^Fli^) (J(X)) dX- 

-^Y.^kj <l''{X)q^x{X) {Fi;,)Fli^) ( J(X)) dX- 



k>0 



fc>0 

]^ f2TT f2TT 

(2'V •'^ -^ fc>0 
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-e / g'^(X)g^(X)— — r.-TY. ^kFi;,) ($(.„)(^ + s(X), f/(X)), . . .) x 



x^Fj)($M(^ + s(X),f/(X)),...j rf'^^rfX + O(e^) (3.14) 

where we used the integration by parts for the generahzed functions. 

We can see now that in the first term of the expression above in the both 
regions X >Y and X < F we have the convolution with respect to X and Y 
of the "slow" functions g3s-(-'^)Q'y(^) with the rapidly oscillating expression 



(i^r.) {^[^n){e + s(X, e), J(X)), . . .) FJ) (<|.(,„)(^ + s{Y, e), J(F)), . . .)) 

in the main order of e. Here (. . .) means the averaging with respect to 
phases 9°'. So as the small AX and AF lead in the main order of e to 
the changes of phase difference equal to -k'^{J{X))^X + 0((AX)^) and 
\k'^{J{Y))/XY + 0{{/XYf) it is not very difficult to see that in the sense 
of "generalized" limit (i.e. in the sense of the convolutions with the "slow" 
functions of X and Y) we can replace this oscillating expression just by its 
mean value 

Y.^^(F^u)){J{x)){Fiy;){j{Y)) 

k>0 
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where (. . .) means the averaging on the space of m-phase solutions. 

As for the last term of (3.14) we remind that its sum with the expression 
arising from the first term of the local part in (3.13) 

I f2n f2n 



dXq^{X)q^{X)j^j^ ...j^ d^-e[ATi^^^n){e + s{X),U{X)),..) 



d 



-J^ekFC,) {^(^n)iO + s(X), f/(X)), . . .) e—F(t,) ($(,„) (^ + s{X), f/(X)), . . . 
is equal according to (3.1) to 

in the leading order of e. 
So we can write now: 

{J q^iX)riX)dX, I q'^iY)J'^iY)dY}\M' = 
[<l''{X)q'xiX){A7) (J(X)) - g5^(X)g^(X) Y: e.(i^wi^5)) {J{X)) + 



k>0 



+q''{X)q^{X)dx {Q^n (AX)) - Ee^(^W^S)) iJiX))]]dX+ 

\ k>0 J 

+e //E efcgx(^)(i^r.)) (AX)) v{X - Y){Fli^) (J(r)) q^iY)dXdY + 0{^) 
•^ "^ fc>0 

After the integration by parts (in the sense of generalized functions) we 
can write this expression in the following "canonical" form: 



\M' 



{J q''{X)r{X)dX, J q^{Y)J^{Y)dY}l 

(^r)(^(X))+Ee.((Fr,)Fr,))(J(X))-(Ff,))(J(X))(Fj))(J(X)))|x 

xg"(X)g^(X)rfX+ 
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fc>0 



which corresponds to the bracket 

{r{x),j^{Y)}\M' = 

= e ({A7){J{X)) + 5: e, ((F('^,)FJ)) - (Ff,))(Fj))) (J(X)) ] 5'(X - Y) + 

\ k>0 J 

,, 2 e.»^.(A- - n»^ . O,^, (3.15, 

fe>0 ^^ ^-^ 

for the functionals J'^{X). 

So according to Lemma 2.4 and the remarks above we obtain for the Dirac 
restriction on A^' 

{e*-{X),J^iY)}n = Oie) (3.16) 

and the relations (3.15) for the brackets {J'^(X), J^{Y)}d in the coordinates 
J(X), ^*(^) and ^o(Xo). 

It is evident also that the Dirac brackets { J'^(X), J^{Y)}j:) on A1' do not 
depend at any order of e on the common initial phase ^o(-^o) because of the 
invariance of J'^(X), the bracket (2.40) and the submanifold A1' with respect 
to the common shifts of 6". 

The dependence of { J'^(X), J^(F)}b on J{X), Jx{X), 9*xiX), ... is reg- 
ular at e ^ and as can be easily seen from (3.15) we have not any depen- 
dence of ^0 ill the first order of e. 

So it is easy to see now that the Jacobi identities for the bracket {...,.. .}d 
on Ai' with coordinates J{X), Oq{X) and ^o(-^o) written for the fields J'^(X), 
J'^(Y) and J^{Z) at the order of e^ coincide with the corresponding Jacobi 
identities for the bracket (3.2). 

So we proved the Jacobi identity for the bracket (3.2). 
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The skew-symmetry of the bracket (3.2) is just the trivial corollary of the 
skew-symmetry of (2.40). 

We now prove the invariance of the bracket (3.2) with respect to the 
choice of the integrals I'^. The proof is just the same as in the local case and 
we shall just reproduce it here. 

Under the condition (IV) (which is the "generic" requirement) we have the 
unique restriction of the Poisson bracket (2.40) on Ai' with the coordinates 
J(X), 6q{X), ^o(-^o) in the form of formal series at e ^ 0. 

So the two restrictions of (2.40) obtained in the coordinates 

and 

irix),enndoiXo)) 

corresponding to the sets {I'^} and {f } (satisfying (A)-(C)) respectively 
should transform one into another after the corresponding transformation of 
coordinates 

rix) = ]io){J{x)) + E ^%)[J^ di]{x) 

k>l 

on M'. 

Now we note that the leading term of (3.15) coinciding with the bracket 
(3.2) transforms according to the transformation 

r{x) = j^,){j{x)) 

which corresponds to the substitution f/'^(X) = U'^{U{X)) on Ai' view the 
relation (2.35). So we obtain the second part of the theorem. 

Theorem 3.1 is proved. 

Remark. 

From Theorem 3.1 it also follows in particular that the procedure (3.2) 
is insensitive to the addition of the total derivatives with respect to x to 
the densities T"^{ip, ^Px, • • •)• This fact however can be also obtained from the 
elementary consideration of the definition of bracket (3.2). 
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and 



Theorem 3.2 

The Hamiltonian functions 

H" = fu''{X)dX 

H = J{VH){U{X))dX 

generate view (3.2) the local commuting flows which give us the Whitham 
equations for the systems (2.11) and (1.6) respectively. 

Proof. 

It is easy to check by the direct substitution that any of H^^ generates the 
"conservative" form 

of Whitham's system for the corresponding flow (2.11). It is easy to see also 
that this flow conserves any of H'^ and so that any H^ and H^^ commute view 
the bracket (3.2). The same property for the Hamiltonian function H (and 
also for the integral of the averaged density of any local integral / commuting 
with H and I'^ and generating the local flow view (1.4)) can be now obtained 
from the invariance of (3.2) with respect to the set {Z*^} since we can use 
the Hamiltonian function H in the form of (1.2) as the one of the integrals 
instead of any of I^ . 

Theorem 3.2 is proved. 
We can also see that the functionals H'^ give us the conservation laws for 
our Whitham system. 

From the Theorem 1.1 it follows also that the flows 

commute with all the local flows generated by local functionals / h{U)dX 
in the Hamiltonian structure (3.2) and it can be also seen that they give us 
Whitham's equations for the corresponding flows (1.7). 
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It can be easily seen also that the described procedure can be applied by 
the same way to the brackets (1.1) written also in the "irreducible" form and 
not only in the "canonical" one. 

The author is grateful to B.A.Dubrovin, S.P.Novikov, V.L.Alekseev, 
O.I.Mokhov, M.V.Pavlov and E.V.Ferapontov for fruitful discussions. 
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